Wetting in electrolyte solutions 



o 

(N 



oo 



o 



I 

C 

o 
o 



> 

00 

in 
o 

rn 
O 
m 



x 

S3 



Ingrid IbagonQ Markus BierQ and S. Dietrich 

Max-Planck-Institut fur Intelligente Systeme, Heisenbergstr. 3, 
70569 Stuttgart, Germany, and IV. Institut fur Theoretische Physik, 
Universitat Stuttgart, Pfaffenwaldring 57, 70569 Stuttgart, Germany 
(Dated: March 13, 2013) 

Wetting of a charged substrate by an electrolyte solution is investigated by means of classical 
density functional theory applied to a lattice model. Within the present model the pure, i.e., salt- 
free solvent, for which all interactions are of the nearest-neighbor type only, exhibits a second-order 
wetting transition for all strengths of the substrate-particle and the particle-particle interactions for 
which the wetting transition temperature is nonzero. The influences of the substrate charge density 
and of the ionic strength on the wetting transition temperature and on the order of the wetting 
transition are studied. If the substrate is neutral, the addition of salt to the solvent changes neither 
the order nor the transition temperature of the wetting transition of the system. If the surface charge 
is nonzero, upon adding salt this continuous wetting transition changes to first-order within the wide 
range of substrate surface charge densities and ionic strengths studied here. As the substrate surface 
charge density is increased, at fixed ionic strength, the wetting transition temperature decreases and 
the prewetting line associated with the first-order wetting transition becomes longer. This decrease 
of the wetting transition temperature upon increasing the surface charge density becomes more 
pronounced by decreasing the ionic strength. 



I. INTRODUCTION 

Wetting transitions are surface phase transitions which 
occur whenever a phase C intrudes at the interface be- 
tween two phases A and B, with either A, B, and C 
in thermodynamic coexistence or with A as a specta- 
tor phase and B and C in thermodynamic coexistence. 
As an example, A is an inert substrate and B and C 
are the gas and the liquid phase, respectively, of a sim- 
ple fluid. The thickness of the intruding liquid film can 
be either finite (incomplete wetting) or macroscopically 
large (complete wetting) upon approaching gas-liquid co- 
existence along an isotherm. The transition at two-phase 
coexistence from incomplete to complete wetting occurs 
at the wetting transition temperature T = T w . It can 
be either continuous (second-order), in which case the 
film thickness diverges smoothly as T — > T w along two- 
phase coexistence, or discontinuous (first-order), imply- 
ing a macroscopically large jump of the film thickness 
from a finite value below T w to a macroscopically large 
one above T w . In the surface phase diagram a first-order 
wetting transition has a prewetting line associated with 
it which is connected tangentially to the gas-liquid coex- 
istence line at T w , extends into the gas phase region, and 
ends at a critical point. The richness of wetting phenom- 
ena has been covered by various reviews [H-Q. 

So far, to a large extent, wetting studies have been de- 
voted to fluids composed of electrically neutral molecules. 
However, for numerous real systems the presence of ions 
is either of crucial importance for wetting phenomena, 
such as electrowetting @, or unavoidable because many 
substrates release ions once they are brought into contact 
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with polar solvents Q ■ For example, electrowetting refers 
to the change of the substrate-fluid intcrfacial tension as 
a response to an applied electrostatic potential difference 
between the substrate and the fluid bulk. This effect 
offers numerous applications in devices based on the ma- 
nipulation of tiny amounts of liquids, such as microti u- 
idic devices H, H[ . Theoretical studies of those systems 
started back in 1938 when Langmuir developed a model 
to determine the equilibrium thickness of water layers 
on planar surfaces in contact with undersaturated water 
vapor, based on the calculation of the repulsive force be- 
tween two plates immersed in electrolyte solutions [Ioj |. 
The typical values for the equilibrium layer thickness as 
predicted by Langmuir's formula were confirmed exper- 
imentally [111 ] and the experimental data were used to 
analyze the effect of various contributions to the dis- 
joining pressure onto the stability of the wetting films 
hj. Some years later Kayser generalized Langmuir's 
model for the equilibrium thickness of wetting layers to 
liquid mixtures of polar and non- pola r components in 
contact with ionizables substrates |13| ; in contact with 
the wetting liquid these substrates donate ions to the 
liquid which act as counterions to the emerging opposite 
charge left on the substrate with overall charge neutral- 
ity. This analysis was followed up by including the effect 
of added salt the ions of which do not stem from the 
substrate [lij]. These papers did not address the issue 
of wetting transitions at coexistence but rather focused 
on the thickness of the wetting layer and the behavior 
of the disjoining pressure. For wetting films of solvents 
without added salt, i.e., with counterions only, Langmuir 
[Toj and Kayser [l3| found that the film thickness / in- 
creases as I ~ (A/i) -1 / 2 , with A/x = fi co — /i, as the 
chemical potential /x approaches its value /i co at coexis- 
tence from the vapor side (p < /x co ). In contrast, wetting 
films without ions and at neutral substrates but with van 
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der Waals interactions (which are not taken into account 
in our model) lead to I ~ (A/z)" 1 / 4 or I ~ (A/i)- 1 ^ 
depending on whether retardation effects are taken into 
account or not, respectively @). In the case that the ef- 
fect of added salt dominates van der Waals interactions 
Kayser [3 found I ~ ln(A/i) as it holds for short-ranged 
interactions. 

Only recently theoretical investigations concerning 
wetting transitions of electrolyte solutions at charged 
solid substrates have emerged fl5l - [l7j . In Ref . [HI the ef- 
fect of adding ions onto the wetting behavior of the pure 
solvent was studied by using Cahn's phenomenological 
theory [Ml for the solvent combined with the Poisson- 
Boltzmann theory for the ions. This model does not take 
into account the solvent particles explicitly, neglecting 
the coupling between solvent par ticles and ions. On the 
other hand, the model in Ref. [16| takes all three types of 
particles (solvent, cations, and anions) explicitly into ac- 
count in terms of hard spheres of different diameters with 
a Yukawa attraction between all pairs and the Coulomb 
interaction between ions. The model was studied by us- 
ing Rosenfeld's density functional theory [l8|, [l9[ com- 
bined with a mean-field approximation for the Yukawa 
and the electrostatic interactions. Within this model, 
the polar nature of the solvent molecules was ignored; it 
was included in a subsequent article by the same authors 
in which the solvent particles were represented by dipolar 
hard spheres [I?) ■ However, for technical reasons, the nu- 
merical analyses of these continuum models in which all 
three types of particles are treated explicitly on a micro- 
scopic level were limited to small system sizes. Therefore 
Refs. [lj| Qj} focused on the case of strong screening 
of the Coulomb interactions which is provided by large 
ionic strengths, i.e., large ion concentrations. However, 
the approaches used in Refs. [l5l - [l~7| are not reliable for 
large ionic strengths due to the use of Poisson-Boltzmann 
theory for the electrostatic interactions which has been 
proved to be valid only for low ionic concentrations [20| . 

In order to overcome these problems we introduce a 
lattice model for an electrolyte exposed to a charged 
substrate which takes into account all three components 
via density functional theory and offers the possibility 
to study significantly broader interfacial regions. In Sec. 
[IT] we introduce this model and the approximate density 
functional. In Sec. IIIII wc present our results for the bulk 
properties and the wetting phenomena for both the salt- 
free solvent and the electrolyte solution. We conclude 
and summarize our main results in Sec. HVl 



II. 



MODEL AND DENSITY FUNCTIONAL 
THEORY 

A. Model 



We study a lattice model for an electrolyte solution 
in contact with a charged wall. The solution consists of 
three components: solvent (0), anions (— ), and cations 



(+). The coordinate perpendicular to the wall is z. The 
region above the wall, accessible to the electrolyte com- 
ponents, is divided into a set of cells the centers of which 
form a simple cubic lattice {r} with lattice constant a. 
The volume a 3 of a cell corresponds roughly to the vol- 
umes of the particles, which are assumed to be of similar 
size. The centers of the molecules in the top layer of the 
substrate form the plane z = 0. At closest approach the 
centers of the solvent molecules and ions are at z = a. 
The plane z = a/2 is taken to be the surface of the planar 
wall. Each cell is either empty or occupied by a single 
particle. This mimics the steric hard core repulsion be- 
tween all particles. Particles at different sites interact 
among each other via an attractive nearest-neighbor in- 
teraction of strength u which is taken to be the same for 
all pairs of particles. In addition, ion pairs interact via 
the Coulomb potential. The solvent particles are taken 
to carry a dipole moment. 

The wall attracts particles only in the first adjacent 
layer via an interaction potential of strength u w which 
is the same for all species. In addition it can carry a 
homogeneous surface charge density a = aea which is 
taken to be localized in the plane z = a/2 and which 
interacts electrostatically with the ions; e > is the el- 
ementary charge. Since we focus on the influence of the 
ions onto wetting phenomena we refrain from consider- 
ing the more realistic, long-ranged van der Waals forces 
which are known to be relevant for wetting transitions 0] • 
Within the mean-field theory we shall use, the choice of 
nearest-neighbor interactions provides a significant com- 
putational bonus which we want to exploit in favor of our 
core concern stated above. 

The corresponding lattice-gas Hamiltonian for this sys- 
tem reads 
1 
2 
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r i 

where rij(r) are occupation number variables, which are 
either or 1 according to whether the cell at the discrete 
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position r = (rii,z > a) = (x,y,z > a) = (ma,na,pa) 
with m,n G Z,_\m\ < M/2 and \n\ < N/2, and 
p = 1,2,3, ••• ,L is empty or occupied by a particle 
(there is no double occupancy); i,j = 0,+,—, eg, is 
the particle charge with q = and q± — ±1; mi(r) is 
the particle dipole moment at r (we consider the typical 
situation of a polar solvent and of ions without perma- 
nent electric dipoles, i.e., m± = 0); w (|r — r'|) = — u for 
nearest neighbors (u > corresponds to attraction) and 
w (|r — r'|) = beyond; —u w is the strength of the at- 
tractive (u w > 0) substrate potential acting on the first 
layer z = a. For the charge density p(r) = aS(z — a/2) 
on a substrate with radial extension Rq the electro- 
static potential is given by (j>(v) = j ^^' i^fezpi = 
£- o WRl + {z- a/2) 2 -\z- a/2\) -> -^z + const, for 
i?o ^> |.z — a/2\ and 2 > a/2. In this regime of be- 
ing close to the charged wall the electric field is uniform 



21| . Therefore the actual position of the charged wall 
enters the electrostatic potential, and thus the Hamilto- 
nian in Eq. ([TJ, only via an irrelevant additive constant. 
The potential energy of a dipole moment nij(r) in the 
electric field E(r) = V^(r) — > T^&z — const, of the 

surface charge is given by — m,(r) • E. In Eq. (p} we 
consider only charge neutral configurations {^(r)}, i.e., 
J2 (n+(r) - n_(r)) = -MNa with r G V = MNLa 3 . 

r 

For weak external electric fields the polarization is ex- 
pected to exhibit a linear response behavior [2l[ . In this 
case, it has been shown that the relative permittivity e 
of microscopic models like the one in Eq. (JTJ) can be 
expressed in terms of molecular properties such as the 
dipole moment and the polarizability [13, . In order 
to simplify our model, the polar nature of the solvent 
is taken into account effectively via the relative permit- 
tivity e of the electrolyte solution which is assumed to 
depend on the solvent configuration no (r) but not on the 
configuration of the ions n± (r) because the orientational 
polarization, i.e., the polarization due to the permanent 
dipoles of the solvent molecules, is the dominant contri- 
bution to the total polarization. In this case Eq. ([T]) 
reduces to (see, c.f., Eqs. (J7J) and (|9])) 



r,r' i,3 



+ 2 J d 6 r*<t>(r*)Q(r*) 



(2) 



where Q(r*) = £ ft<0*) + o5{z* - a/2) is the local 

i 

charge density where n*(r*) = n^(r) for all r* <G (aR) 3 
and r g (aZ) 3 with max (\x* — x\, \y* - y\, \z* — z\) < 
a/2; 4>(r*) is the electrostatic potential which can be ob- 



tained by solving the Poisson equation 

- £ „V-[ £ K(r*))V0(r*)] = Q(r*, [n* ± ]), r* G (aR) 3 nV, 

(3) 

where V is the volume of the fluid. For general permit- 
tivity profiles £(rig(r*)) no closed solution <j)(r*) of Eq. 
(|3|) as a functional of e(ng(r*)) and Q(r*) is known, i.e., 
for each configuration {rii(r)} the evaluation of Eq. 
requires to solve the differential equation ([3]) anew. It 
has been proven, that models including charges as in Eq. 
(|2"| possess a proper thermodynamic limit for sequences 
of finite-sized systems, which is independent of the shape 
of the container, provided that globally charge neutral 
configurations {rij(r)} are considered [2J, H^. Since the 
thermodynamic limit is performed for sequences of finite- 
sized systems the electrostatic potential <j>(r*) in Eq. ([3]) 
vanishes at infinity (|r*| — > oo) plj . 



B. Density functional 

With a given expression for e(no(r)) (see, c.f., Eq. 
(fl5]l). Eq. ([2]) can be used directly for numerical anal- 
yses such as Monte Carlo simulations, provided an ef- 
ficient method to determine the electrostatic potential 
cf>(r*) for arbitrary permittivity profiles e(no(r)) becomes 
available (see for example Ref. [26[ for recent efforts in 
this direction). We leave this challenging task for future 
studies. Here, we consider a suitable mean field approxi- 
mation which can be formulated as to minimize a grand 
canonical density functional fi[{pi(r)}] [27j of continuous 
and dimensionlcss occupation number distributions pi (r) 
such that at the minimum pi(r) = /0^ 9 (r) approximates 
the thermal average (rij(r)). 

Application of the Bragg- Williams Approximation [2E 
[3T| to the model Hamiltonian in Eq. leads to the 
following grand canonical density functional: 



/?n [{*(*)}] = £ 



^Pi(r) In pi (r) 

i 

l-^ Pi (r))ln(l-^ ft (r)) 

i 3 



r,f' i,j 



2nl R / d A f 



_ (D 



r 1 

2 



v e{pl{v*)) 



(4) 



where /3 = (fc B r) _1 is the inverse thermal energy 
and pi is the chemical potential of species i, Ib 
/ ;> a = e 2 /3/(47reo) is the Bjcrrum length in vacuum, 
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r = r/a are the dimensionless lattice positions, r* = 
r*/a, p*(f*) = pi(f) for all r* e K 3 and r e I? 
with max (\x* - x\, \y* - y\, \z* - z\) < 1/2. The ac- 
tual number densities of the components are given by 
p%(f) = pi(r)a~ 3 . Charge neutrality demands J2iP+(?)~ 



Aa 2 



MNa is 



p-(r)] = -Act where A = MN 
the substrate area and a = <r/(ea~ 2 ); this constraint 
is implemented via a boundary condition for D (see, c.f., 
Eq. (|T4|) ). The first two terms of Eq. ((4]) represent the 
ideal gas or entropic contribution F it i to the Hclmholtz 
free energy functional F[{pi(f)}} = F id [{pi{r)}\ + 
F e x [{Pi(r)}]; the third and the fourth term represent 
the non-electrostatic contribution to F ex [{/Oj(r)}], which 
follows from the first and second term in Eq. and 
turns out to be equal to the random phase approxima- 
tion (RPA) within density functional theory [23]. This 
approximation is justified because it has turned out that 
RPA is reliable in the present situation of vanishing con- 
trast between the non-electrostatic interactions of the 
three species (Hf. The last term is the electrostatic en- 
ergy Using SI units, the electrostatic field energy density, 
which enters into Eq. (0}, is given by plj 



-ED 

2 



1 D 



ID 2 - 2 



= 2 7 rfc i3 TZB^3 (5) 
2 e e 2 e sa^ ea A 



and 



1e D = - l -V4> ■ D = \4> (V • 6) - iv • (0D) , (6) 

where E = — = — is the actual electric field, 6 is the 

electrostatic potential and D = Dea -2 is the actual elec- 
tric displacement generated by the ions and the surface 
charge density a = aea~ 2 , satisfying Gauss' law [2l[ 



V-D = Q(r*), 



(7) 



so that (V = aV) 

V • D (r*, [p* ± ]) = £>p*(r*) + - 1/2). (8) 



Due to Eq. ((HJ), the electrostatic contribution to the 
functional can be written as 



F e , = I I d A i 



r*) V-D 



)D 



(9) 



where the last term leads to a vanishing surface contribu- 
tion [2lj . because the thermodynamic limit is performed 
for sequences of finite-sized systems. Using Eq. ([7]) ren- 
ders the last term in Eq. ([2]). 

Because the substrate potential depends only on z, the 
minimum of /3Q [{pi(f)}] lies in the subspace of distribu- 
tions Pi(j) which depend on z only. Therefore we write 
Eq. @ for the special case pi(f) = Pi(z), i.e., 



m [{Pi(z)}] 



2=1 I i 

+ (l-J2pi(z)) In (l 

i J 

~ puJ2{pi(z)pj(z + !) + 2pi{z) Pj (z)) 

ij 

-fiu w ^2pi(z)Sz,i - P^fMpijz) > 

i i ) 

L+1 />_*(D(z*,[ p * ± ])) 2 



2nl B / dz* 

h/2 e{p* (z)) 



(10) 



where A = Aa 2 is the substrate area so that AL is the 
volume of the fluid (L = La), and pi(L + 1) = 0. 
Gauss' law (Eq. ((8])) reduces to 

dD(z* > 1/2, [p* ± ]) ^ 

i 

(ii) 

where the last term in Eq. © appears as a boundary 
condition to Eq. (fTT]) : 



D(z* = l/2,\p* ± ])=a. 



(12) 



Since p* ± E [0,1] arc bounded, i.e., the densities p± do 
not exhibit <5-likc singularities, the boundary condition is 
determined entirely by the surface charge. 

The density profiles p± (z) have to fulfill global charge 
neutrality, i.e., 



^2\p+(z)-p-(z)]+a = 0, 



(13) 



2=1 



which according to the integrated Eq. (|lll) is equivalent 
to 



D (z* = L + 1/2, [p* ± ]) = 0. 



(14) 



The relative permittivity s(z*) is taken to depend lo- 
cally on the solvent density Pq(z*) through the Clausius- 
Mossotti expression [2l| 



<Po(z*)) 



1 



2a 
3£o 



Po(z*) 



(15) 



where a is an effective polarizability of the solvent 
molecules. In the following its value is chosen such that 
e = 60 for po = 1; this choice corresponds to a mean 
value for liquid water along the liquid- vapor coexistence 
curve. 

As for a lattice model Eqs. ([2]) and (fT0|) do not include 
the kinetic energy. The latter requires an off-lattice de- 
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scription which leads to a density independent contribu- 
tion to the chemical potential of species i so that 

Vi.phys = k B Tln(piAf) + p eXl 

= k B T\n{^A*) + p ex , 

a- 3 (16) 

= k B Tln(p t ) + p ex + 3fc B Tln(A l /a), 
= fii + 3fc s Tln(Ai/a), 



where A, = hj \J2ixmik B T is the thermal wavelength , mi 
is the particle mass, and p ex is the excess chemical poten- 



tial over the ideal gas contribution. This gives rise to a 
density independent difference between the actual physi- 
cal chemical potential pi^ p h ys and the chemical potential 
Hi of the lattice-gas model: p-i, p hys — Ma = 3k B Thi(Ai/a). 

C. Euler-Lagrange equations 

In order to obtain the equilibrium configuration, the 
density functional in Eq. (|10[) has to be minimized under 
the constraints given by Eq. <jT2]) and Eq. (TT3J The 
variation of Eq. (flQ|) reads: 



m [{ Pi (z)}] 



\npi(z) - (3p, - In ( 1 - E^( 



E EM* 



/37^(^(z)^(z + l) + ^(2)^(2+1) + 2Spi(z) Pj (z) + 2p i (f)5p i (^ 



2tt// 



Airh 



+1 '\_ ^D{z\[p* ± ])) 2 . ..... 

dz VTT7^^2~ e WW ))2^**.o^i(* ) 
5 -i/2 (e(pg(**))) T 

^* ^/V^tP SD(z*) - K V 5 s ,i5pi(z) 

2-1/2 e lPol z 77 , 



lnpj(z) - Pm - In (l ~Eft( 



2=1 I j 



1j 

-j8«E + E Pi(z') 6 *,z'+i + 2 Pj(z) + 2pj(«)) - ^«w*s,i 



-27rZ/ 



i+1 / 2 d _.(^(^,K])) 2 ^^. 



5-1/2 



(e(p8(**)))' 



£ '(pS(r))^, <5p*(r) 



2+1/2 



-1/2 



dz*(/>'(z*)8D(z*) 



(17) 



I 

where </>(z*) = f3ecf>(z*) is the dimcnsionless electro- zG {1,...,L} 
static potential which fulfills 



D(z*) = e eE(z*) = -e Q e 



1 d 



d<f>(z*) 
dz* ' 



ea z D{z*) = -e e- — 



D(z*) 



4ttZ, 



a dz* \ fie 

>'(r). 



(18) 



hipi(z) - p* - Puy,8z,\ - hi ( l - Eft( 



2+1/2 



2-1/2 



dr«^(r) 



2tt// 



*-+v^(^ M/ 4])) 2 ^ 



2-1/2 



(eOtfC**)))' 



e' (p5(z*))5i,o = 



(19) 



Upon integrating by parts the last term in Eq. (|17l) . 
by using Eq. (TT5J so that 8D(z* = 1/2) = and Eq. 
Cg so that (5L>(z* = L + 1/2) = 0, with &D'(S*) = 

duc to E q- dHJ), and (5p*(z*) = 8p t (z) for with i,j = 0,+,-, where q l e is the electric charge of 

all z* G K and leZ with max(|z* - z| < 1/2) we obtain component i and T* — ^ is the reduced temperature 

the following three coupled Euler-Lagrange equations for and p* = f3pi. At the wall the convention pj(0) = is 
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used. The integrals in Eq. (fT9|) are approximated by 

i-z+l/2 

/ dz*f(z*) « 
Jz-1/2 

((z + 1/2) - (z - V 2 )) / I g J = 

= /(*). (20) 

For given chemical potentials pi these coupled equa- 
tions can be solved numerically by an iterative algorithm. 
The values of the chemical potentials pi considered here 
correspond to those for the bulk gas phase of the system. 
For each iteration the electrostatic potential 4>(z*) must 
be calculated by solving Poisson's equation (see Eqs. (fTTI) 
and QD)) 

A.(e( p *(z*))^(z*)) = -Airlz^UpKn (21) 

i 

ensuring global charge neutrality at each step. 

D. Wetting films 



of 647 K for water. In our units 1 mM = 10~ 3 mol/L 
corresponds to pj = pid 3 = 3.9 x 10 -5 . 

For our calculation we have used values for the reduced 
surface charge density a in the range between and 10~ 2 . 
For a = 4A the latter value corresponds to 1 pC/cm 2 . 
Such values are within the range of measured surface 
charge densities of silicon nitride at two different con- 
centrations of the background electrolyte NaCl (1 mM, 
10 mM) determined by potentiometric pH titration [33j . 
which is a common method to determine the unknown 
concentration of an identified substance and to estimate 
the surface charge of a solid by comparing the titration of 
the solution with solid against the titration of the same 
solution without solid. 

These consideration indicate that the values of the 
reduced substrate surface charge densities a and ionic 
strengths / considered in the following are within the 
range of values for which Poisson-Boltzmann theory, i.e., 
mean-field theory for the electrostatic interaction, shows 
quantitative agreement with corresponding Monte Carlo 
simulations [20]. The former is essentially identical to 
the theory used to describe the ions in Eq. (|TU|) if one 
neglects the effect of nonzero ion size, which is weak for 
the considered dilute electrolyte solutions. 



The wetting behavior can be transparently inferred 
from the constrained surface contribution f2 s (I) := 

}] — to the grand potential Q, where is 

the bulk contribution to the grand potential and the den- 
sity profiles pf' are the solutions of the Euler-Lagrange 



equations (|T 
fined as 



for a prescribed film thickness I = la de- 



/ = 



P0,1 - PO.g 



(22) 



where T = J Q °° dz (po{z) — /5 (oo)) = Ta~ 2 is the excess 
adsorption (or coverage) of the substrate by the solvent 
and po,i and po, g are the corresponding bulk number den- 
sities of the liquid and the gas phase, respectively. In 
order to obtain fl s (I) by using a Lagrange multiplier we 
have minimized [pi(z)] under the constraint 

00 

^2(po(z) - po.b) = r = l{po,\ - po, g ), (23) 



where po,fc is the number density of the bulk gas phase in 
units of a -3 . 



E. Choice of parameters 

If one chooses the lattice constant a to be equal to 4A, 
the maximal density 1/a 3 lies between the densities for 
liquid water at the triple point and at the critical point. 
Accordingly, the choice Ib = 400 corresponds to T w 
417 K. This temperature lies between the triple point 
temperature of 273 K and the critical point temperature 



III. RESULTS AND DISCUSSION 

A. Bulk Phase Diagram 

In the bulk, the number densities pi of the fluid are spa- 
tially constant and from the requirement of local charge 
neutrality it follows that p + = p_ = /, where I is the 
so-called ionic strength for monovalent ions. Under these 
conditions the density functional given by Eq. (fTUf re- 
duces to 



V 



po(lnp - Mo) + J(21nl - p}) 



(l-po- 21) In (1 - po - 2/) - — (po + 2l f 



(24) 



where p} = p* + + p*_ and V — V/a 3 (V is the volume of 
the fluid). The last term in Eq. (TITj|) vanishes because 
in the bulk D = due to Eq. (TTTj) . The Euler Lagrange 
equations (|19p read 



In po - Po - In (1 - Po - 21) - — ( Po + 21) = 



2lnl -p*j -21n(l 



Po -2I)-—(po + 2I) = 0. 



(25) 



For a given ionic strength I = p± in the liquid phase 
of the solution, the liquid-gas coexistence curves, i.e., the 
solvent density in the liquid phase of the solution and the 
coexisting densities of the ions and of the solvent in the 
gas phase of the solution, are determined by the equality 
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-1 




FIG. 1. Bulk phase diagram no, co (T) of liquid-gas 

coexistence according to Eq. (|24|) in the fig — T* plane for 
the salt-free (J = 0) case of a pure solvent. If the wetting 
transition temperature T£ is above the triple point T t * ~ 0.21 
(for water), three types of paths (A), (B), and (C) are used 
to study the wetting behavior of our model. (A) is a path 
along gas-liquid coexistence on the gas side whereas along 
the paths (B) and (C) two-phase coexistence is approached 
along isotherms leading to incomplete (B) and complete (C) 
wetting, respectively. 

of the chemical potentials po and pi and of the pressure 
p: 

Mo[{p J (9) },n = Mo[{p J (1) },n, 

Mi[{4 s) },T*]= Mj [{p«},T*], (26) 
p[{p ( ? ) },n=p[{p ( P},T*}. 

For I~0 the resulting phase diagram can be deter- 
mined analytically and is plotted in Fig. [T] The reduced 
critical temperature is T*(I = 0) = 0.5 and the criti- 
cal number density is po lC (I = 0) = 0.5. For 1^0 the 
binodal curves are determined numerically and the criti- 
cal points are obtained by determining the maximum of 
the corresponding spinodal curves. Within the present 
model the reduced critical temperature T* is indepen- 
dent of I whereas po. c (I) = 0.5 — 21. In agreement with 
experimental evidence [341 ] the shift of the binodal curves 
is negligibly small for ionic strengths up to 10 mM, i.e., 
I < 3.9 x HT 4 . 



B. Wetting 

1. Salt-free solvent 

We first consider the case I = 0, in which our model 
reduces to the lattice-gas model studied by Pandit et al. 
[35L HH . In that case, the Euler-Lagrange equations in 



Eq. ([19]) reduce to 

lnpo(-z) - hi[l - po(z)] - Mo _ P u wSi,z 

-3^[4p(z) + p(« + l)+/>(z-l)] = 0, (27) 

and the ratio u w /u = 3T*f3u w controls the wetting and 
drying transitions. For u w /u > 1 the substrate is so 
strong that it is already wet at T* = 0; in the range 
0.5 < u w /u < 1 there is a wetting transition at > 0; 
and in the parameter range < u w /u < 0.5 a drying 
transition occurs. Depending on the value of the ratio 
u w /u one observes layering transitions, i.e., one can dis- 
tinguish the number of discrete layers which arc forming 
upon reaching thick films. The transition from n to n+ 1 
layers is first order and shows up as a jump in the film 
thickness I. The loci of these discontinuities are layer- 
ing transition lines, each ending at a critical point T* n . 
For large n, T* n approaches the roughening transition. 
However, within the present mean-field theory T* n ap- 
proaches T* . Since layering transitions should only oc- 
cur along or near the melting curve or the sublimation 
line, these layering transitions are a special feature of 
the lattice-gas model used to describe the liquid and gas 
phases Q. 

We have carried out calculations in the parameter 
range 0.5 < u w /u < 1. A wider range of the pa- 
rameter u w /u was studied thoroughly by Pandit et al. 
[35l . . Figure [2] shows the effective interface potential 
to(l) = £l s (l) — J g ,i — 7/, s for three different temperatures 
along a path at coexistence [path (A) in Fig. [T] for the 
rather arbitrarily chosen values u w /u = 0.81 and u w /u = 
0.69. Here j g j and 7; ;S are the gas-liquid and liquid- 
substrate interfacial tensions, respectively, such that by 
construction at two-phase coexistence u)(l — > 00) = 0. 
The equilibrium thickness of the liquid film is given by 
the position of the global minimum of uj(1). If I = 00 
is the global minimum of i7 s (Z) the system is wet. In 
this case, the gas-substrate surface tension is given by 

7g, S = tt S {l = OO) = JgJ + 7;, S . 

In the two cases which we have considered in Fig. 
[21 to(l) exhibits only a single minimum, the position of 
which diverges continuously or via steps of finite size 
as T* — > T*. For T* > T* the position of the min- 
imum is / = 00 and the system is wet. The wetting 
transition is second order and occurs at the temperature 
T* ~ 0.856T* for u w /u = 0.81 and at T* ~ 0.95T* for 
u w /u = 0.69. Within the present model, in which all 
interactions are of the nearest-neighbor type only for the 
pure solvent, the system exhibits a second-order wetting 
transition in the entire parameter range 0.5 < u w /u < 1. 
This observation is compatible with corresponding Monte 
Carlo simulations of the Ising model on a cubic lattice 
[37l . l3a |. However, the order of wetting transitions de- 
pends sensitively on the range of interactions as well as 
on whether a continuous or a lattice model is considered. 
For a continuous analogue of the present model, Pandit 
et al. (36T ] found a second-order wetting transition only 
for 0.5 < u w /u < 0.7 but a first-order one for u w /u > 0.7. 
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proachcd along four isotherms from the gas phase [paths 
of type (B) and (C) in Fig.Q], is plotted in Fig. |3 In 
the case u w /u = 0.81 (Fig. 3(a)) the isotherms exhibit 
vertical steps at the aforementioned layering transitions. 
Above i.e., when the substrate is completely wet at 
coexistence, the isotherms exhibit an unlimited number 
of such steps as Po co (T*) — /ig approaches zero, while 
for T* < there is only a finite number of steps. For 
u w /u = 0.69 (Fig. 3(b)) layering transitions do not oc- 
cur and the film thickness diverges logarithmically for 
T* > T*, while for T* < T* it reaches a finite value at 
coexistence. 



2. Electrolyte solution 
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FIG. 2. Effective interface potential lu(1) — Q s (t) — 7 9 ,z — 7z, s 
at two-phase coexistence as a function of the thickness / = la 
of the adsorbed liquid film for three temperatures in the salt- 
free case (/ = 0) for u„/u = 0.81 (a) and u w /u = 0.69 (b) . In 
both cases u)(l) exhibits only a single minimum, the position 
of which diverges continuously as T* — s- TZ- Accordingly, 
the system undergoes critical wetting at T£ ~ 0.856T C * for 
u w /u = 0.81 and at TZ ~ 0.95T* for u w /u = 0.69. 



Moreover, lattice-gas models with short-ranged 
particle-particle interactions and long-ranged substrate 
potentials were studied by de Oliveira and Griffiths [39[ 
and Ebner [13, S3 . In Ref. HU complete wetting in 
a system with T w = was studied within mean field 
theory. Ebner reported T w = or a first-order wetting 
transition depending on the strength of the substrate 
potential [4(| and studied the same interaction poten- 
tials as the ones used in Refs. [39|, |4(j applying Monte 
Carlo simulations [4l[ . Finally, in systems in which both 
the particle-particle interactions and the substrate po- 
tentials are long-ranged, critical (i.e., second-order) and 
first-order wetting can occur for suitable choices of the 
interaction potentials [42], |43[ . 

The film thickness / = l/a as function of /Zq co (T*)—/j,q, 
when bulk coexistence Hq. co (T*) (see Fig. [I} is ap- 



Within the above concepts we now focus on the influ- 
ence of the ionic strength / = Ia~ 3 and of the surface 
charge density a = <rea~ 2 on the wetting behavior of sys- 
tems with u w /u = 0.81 or u w /u = 0.69. If the substrate 
is neutral (a = 0), the addition of salt changes neither 
the order nor the transition temperature of the wetting 
transition, i.e., there is a second-order wetting transi- 
tion at the wetting temperature as discussed in the 
previous Subsubscc. IIII B II This is expected because 
within our model all particles have the same size, the 
ions have the same absolute charge, and the strength of 
the particle-particle and of the substrate-particle nearest- 
neighbor interactions are the same for all three species. 
Hence local charge neutrality (p+(z) = p~{z)) holds 
due to the exchange symmetry with respect to the ionic 
components. This implies that there is no electric field 
(D(z) = 0). If the surface charge becomes non-zero, the 
order of the wetting transition changes from second or- 
der (er = 0) to first order (<r ^ 0) for all values of the 
charge density a and ionic strength I studied here, with 
a = 2 x 10~ 5 (i.e., a w 0.002 ^C/cm 2 ) as the smallest 
non-zero value considered. This result is in agreement 
with previous studies. The influence of ionic solutes on 
the order of the wetting transition was studied in Ref. 



15 1 by using Cahn's phcnomcnological theory and in Ref. 

16| by using density functional theory for an explicit sol- 
vent model for an ionic solution. Both studies suggest 
that electrostatic interactions favor first-order wetting. 

Figure 0] shows examples of the effective interface po- 
tential uj{1) in the case of non-zero surface charge densi- 
ties, a = 2 x 10 -3 and a = 2 x 10 -4 , for two temperatures 
and at bulk coexistence [sec path (A) in Fig. [Tj. In both 
cases, ui(l) has two local minima. For T* < the global 
minimum corresponds to a thin film whereas for T* > 
the film is macroscopically thick. At the wetting transi- 
tion temperature the two minima correspond to the 
same value of the effective interface potential oj(1). Ac- 
cordingly, at the film thickness jumps discontinue usly 
from a finite value below to a macroscopic one above 
so that the system undergoes a first-order wetting 
transition. If a is decreased the height of the barrier in 
ui(l) at the wetting temperature decreases and the 
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FIG. 3. Film thickness I = l/a in units of the lattice con- 
stant a as a function of undersaturation Ho }CO (T*) — fio for 
the salt-free case (I — 0). Gas-liquid coexistence Ho, co (T*) 
is approached from the gas phase, (a) u v /u = 0.81: for 
T* < = 0.856T* the system is partially wet and, if 
at all, there is a finite number of layering transitions; for 
T* > the isotherms exhibit an unlimited number of lay- 
ering transitions as fjZ,co(T*) — fio — > and the first few lay- 
ering transitions are rounded because for this temperature 
T* > T c * n . (b) Uw/u — 0.69: the film thickness diverges 
logarithmically for T* > T* = 0.95T C *, while it reaches a 
finite value at coexistence for T* < In (b) there are 

no layering transitions. Note that with T*(I) = | one has 

2/3 |"w),co(T*)-Mol 



Mo,co(r*)-Mo 



T* /T* 



minimum close to the wall is shifted to larger thicknesses 
(Fig. IHb)). In the case a = 0, uj(1) has only a single 
minimum, like in the salt- free case (see Fig. [2|, corre- 
sponding to a second-order wetting transition. 

In Fig. [5] the wetting transition temperature is plotted 
as function of the surface charge density for two values 
of the ionic strength and for u w /u ~ 0.81. As a ~ aa 2 /e 
is increased, the wetting transition temperature de- 
creases due to the strengthening of the substrate-fluid 



FIG. 4. Effective interface potential uj{1) at gas-liquid coexis- 
tence as function of the thickness I = l/a of the liquid film for 
u„/u = 0.69, / = 3.9 x (/ = ImM), and a = 2 x 10~ 3 

(a = 0.2/iC/cm 2 ) in (a) and a = 2 x 10~ 4 (a = 0.02/iC/cm 2 ) 
in (b) for two temperatures in each case. The effective inter- 
face potential u)(l) has two local minima (one at I < oo and 
one at I = oo) which have the same depth at T^. Accordingly, 
for both surface charge densities a the system undergoes a 
first-order wetting transition. 



attraction as the substrate is charged up. For a ^ the 
system with a smaller ionic strength / has always the 
lower wetting transition temperature because in this 
case the screening of the electrostatic forces of the sub- 
strate is reduced making them effectively stronger which 
favors wetting. As already mentioned above, within our 
model for a = the wetting transition temperature is 
independent of the ionic strength / = Ia~ 3 . The trend 
is the same for u w /u = 0.69. In the case of first-order 
wetting transitions these results are in qualitative agree- 
ment with Ref. [l6||. However, the off-lattice model used 
therein exhibits also second-order wetting transitions (see 
the discussion above in Subsubscc. IIII B for which T* 
is a non-monotonic function of a. 
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FIG. 5. Wetting transition temperature as a function of 
the substrate surface charge density a — era 2 /e for u w /it = 
0.81. The two types of symbols correspond to distinct values 
of the ionic strength / = la 3 in the bulk liquid phase (• for I = 
3.9 x 10" 5 (/ = ImM) and ■ for I = 3.9 x 1(T 4 (J = 10mM)). 
Filled symbols correspond to first-order wetting transitions, 
while the empty one at a — corresponds to a second-order 
wetting transition, with the corresponding wetting transition 
temperature being independent of I. 



Since the wetting transitions for a ^ are first or- 
der, there is a prcwetting line associated with them. The 
prewetting line is attached tangentially to the gas-liquid 
coexistence line at the wetting temperature T£ and bends 
away from coexistence, marking the loci of a finite discon- 
tinuity in film thickness I = l/a. The discontinuity upon 
crossing the prewetting line gets smaller as one moves 
further away from coexistence and it vanishes at the 
prewetting critical point. Figure [S] shows the film thick- 
ness I = l/a for four different isotherms as a function of 
undersaturation Pq, co (T*) — Hq for u w /u = 0.81 and a = 
2 x 10 -3 (a — 0.2/iC/cm 2 ). The film thickness increases 
for small undersaturation as I ~ ln(/iQ co (T*) — /ij). Ac- 
cordingly, w(7) ~ exp(— 2kI), where k = y / 8nlBl/e(p 1 ) 
is the inverse Debye length (see inset of Fig. [6]). This 
is in agreement with Refs. [14j and 15j for wetting of 
solvents with added salt. In contrast, for wetting films 
of solvents without addition of salt, i.e., with counte- 
rions only, one has I ~ (/i co — /i) -1 / 2 and cj(1) ~ I" 1 
[Tol . [12, El HI ■ In order to obtain this result, Eqs. ^ 
and l|10p have to be modified to consider only solvent 
particles and countcrions but leaving out coions. In ad- 
dition to the finite thin-thick jumps in film thickness I 
when crossing the prewetting line we observe first-order 
layering transitions similar to those found in the salt-free 
case for u w /u = 0.81 (see Fig. [3j). The addition of the 
electrostatic interaction leads to a series of triple points 
where the layering transition lines meet the prcwetting 
line, as shown in the surface phase diagram in Fig. [7J A 
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FIG. 6. The film thickness I = la (Eq. ([22)1 as a function of 
undersaturation ^o, co {T*) — fJ,o along four different isotherms 
for u v /u = 0.81, I = 3.9 x 10~ 5 (i = ImM), and a = 2 x 10~ 3 
(a — 0.2/iC/cm 2 ) exhibits a large but finite jump (corre- 
sponding to more than one monolayer) when the prewet- 
ting line is crossed and small jumps when the various layer- 
ing transition lines are crossed. The film thickness increases 
for small undersaturation as I ~ ln(/^o. co{T*) — fio) where 
M(*, co (r*)-M5 = T^7% [ TO '°° ( ^ ) ~ M0 ]. The inset displays 
the corresponding asymptotic behavior of the effective inter- 



face potential uj(1) ~ exp(— 2kI) where k = \J8ttIb I/(£(Po^ 
is the inverse Debye length. 



similar phase diagram was found by Ebner [40( using a 
lattice-gas model for a one-component fluid in which the 
fluid particles interact among each other via a Lennard- 
Jones (6-12) potential and a fluid particle interacts with 
the substrate via a (9-3) potential. This is also in line 
with the prediction by Pandit et al. [13 for a substrate 
of intermediate strength, i.e., for 0.5 < u w /u < 1, with 
interactions ranging beyond nearest neighbors. 

In the case u w /u = 0.81 and for fixed ionic strength / 
we have studied the prewetting lines for various values of 
the surface charge density a. Figure [5] shows the prewet- 
ting lines for ionic strength / = 3.9 x 10~ 5 (/ = ImM) 
and for four values of a. One can see clearly that as 
a decreases, the wetting temperature rises and the 
prcwetting line becomes shorter. This is in agreement 
with the fact that in the limit a — > the wetting tran- 
sition turns second order. The values of the prcwetting 
critical points for the lines shown in Fig. [S] are given in 
Table [0 



IV. CONCLUSIONS AND SUMMARY 

We have investigated wetting of a charged substrate 
by an electrolyte solution with a focus on the influence 
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FIG. 7. Surface phase diagram for u v /u = 0.81 and a — 
2 x 1CT 3 (a = 0.2^C/cm 2 ) . The full line is the prewetting 
line attached to T* — 0.864 T* (a) and ending at the prewet- 
ting critical point (•). The dashed lines correspond to layer- 
ing transition lines. They end at layering critical points T* n 
(located at the end of the dashed lines without being indi- 
cated separately), which within the present mean- field theory 
accumulate for n — > oo at T* instead of at the roughening 
transition temperature of the gas-liquid interface on the lat- 
tice. 
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FIG. 8. Prewetting lines for four values of the surface charge 
density a = era 2 /e with ionic strength I = 3.9 x 10~ 5 (/ = 
ImM) in the bulk liquid phase and for u w /u = 0.81. The 
locations of the wetting transitions (a) and of the prewetting 
critical points (•) are given in Table [I] 



of the ionic strength I and of the substrate surface charge 
cr on the wetting behavior. First, we have investigated a 
lattice-gas model for the salt- free, i.e., pure solvent (Fig. 
[l} providing a reference system relative to which the in- 
fluence of the electrostatic interaction can be compared. 
The results for the salt-free case are in good agreement 
with previous studies [3]| HH . We have calculated the ef- 
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TABLE I. Prewetting critical points {T* w c ,fiQ pvr c ) for the 
prewetting lines shown in Fig. [8] The ionic strength in the 
liquid phase is I = 3.9 x 10~ 5 (/ = ImM). T„ is the transition 
temperature for first-order wetting. Note that /Jo, C o(Ipw,c) — 

* 2/3 | >0,eo(Tp We )-M0,pw,c ] 

MO.pw.c - T p * w C /T« [ u J • 



fective interface potential ui(l) which facilitates the trans- 
parent identification of the order of the wetting transi- 
tion (see Fig. [2]). Depending on the value of the ratio 
u w /u of the strengths of the substrate potential and of 
the particle-particle interaction, the model can exhibit 
layering transitions when gas-liquid coexistence is ap- 
proached along an isotherm (see Fig. [3]). In the next 
step we have analyzed quantitatively the effects of the 
ionic strength and of the surface charge density on the 
order and on the transition temperature of the wetting 
transition. Concerning the order of the transition we 
have found that electrostatic forces induce a first-order 
wetting transition, even for very small surface charges 
(see Fig. [4}. Within our model, for cr = the transition 
is second order and occurs at the same temperature as in 
the salt-free case. For a fixed ionic strength, the wetting 
temperature decreases with increasing surface charge 
density of the substrate. This is due to the increasing 
substrate-fluid attraction as the substrate surface charge 
is increased. If systems, which differ only with respect to 
the ionic strength /, are compared, the one with smaller 
/ has the lower wetting transition temperature (see 
Fig. [5]). When bulk coexistence is approached along an 
isotherm, in the case of a first-order wetting transition, 
i.e., if cr 0, the model exhibits first-order layering tran- 
sitions in addition to prewetting (Fig. [6]). This leads to 
a series of triple points in the surface phase diagram (see 
Fig. [7]). We have also studied the influence of the surface 
charge density on the prewetting lines. We have found 
that the prewetting line becomes shorter as the surface 
charge density is decreased (Fig. [5]). 

Although our lattice model differs significantly from 
the continuum models used in Refs. (la . Il6j . we have 
arrived at similar conclusions concerning the trend that 
adding ions promotes first-order wetting transitions. Ac- 
cordingly, this result can be considered to be robust. 
Within our approach one is able to study wide intcrfa- 
cial regions and therefore small ionic strengths which was 
not possible within the model studied in Ref. [16j. Our 
study also includes a discussion of prewetting, providing 
a more complete description of the wetting properties of 
electrolytes. In agreement with Refs. [3 HH the growth 
law of the film thickness for complete wetting along an 
isotherm is not changed by adding ions to the solvent, in 
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spite of their long-ranged Coulombic interaction (Figs. which are donated by the substrate and the charge of 
[3] and [5]) . However, if only counterions are considered, which is opposite to that of the wall, the film thickness 

varies as I ~ (pi co - fi)- 1 / 2 [M EE El . 
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